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Abstract 



We consider the nonlinear Dirac equation in one dimension, also known as the Soler model in (1+1) 
m ' dimensions, or the massive Gross-Neveu model: 



idtip^-iad x ip + m/3ip-f(i)*0i))f3ip, tp(x,t)eC 2 , xeR; geC°°(R), m > 0, 



where a, ft are 2x2 hermitian matrices which satisfy a 2 — f3 2 — 1, a/3 + /3a — 0. We study the spectral 
stability of solitary wave solutions 4>u(x)e~ luJt . More precisely, we study the presence of point eigenvalues 
in the spectra of linearizations at solitary waves of arbitrarily small amplitude, in the limit ui — >■ m. We 
prove that if f(s) = s k + 0(s ), k € N, with k > 3, then one positive and one negative eigenvalue are 
present in the spectrum of linearizations at all solitary waves with uj sufficiently close to w = m. This 
shows that all solitary waves of sufficiently small amplitude are linearly unstable. The approach is based 
on applying the Rayleigh-Schrodinger perturbation theory to the nonrelativistic limit of the equation. 

The results are in formal agreement with the Vakhitov-Kolokolov stability criterion. 

Let us mention a similar independent result IGG1 on linear instability for the nonlinear Dirac equation 
in three dimensions, with cubic nonlinearity (this result is also in formal agreement with the Vakhitov- 
Kolokolov stability criterion). 



1 Introduction 

A natural simplification of the Dirac-Maxwell system [Gro66[ is the nonlinear Dirac equation, such as the 
massive Thirring model [Thi5& with vector-vector self-interaction and the Soler model [Sol70 with scalar- 
scalar self-interaction (known in dimension n = 1 as the massive Gross-Neveu model |GN74[|LG75j ). There 
was an enormous body of research devoted to the nonlinear Dirac equation, which we can not cover compre- 
hensively in this short note. The existence of standing waves in the nonlinear Dirac equation was studied in 
|Sol70j , |CV86] , |Mer88] , and [ES95 . Numerical confirmation of spectral stability of solitary waves of small 
amplitude is contained in [Chu07] . The overview of the well-posedness of the nonlinear Dirac equation in ID 
is contained in [PellOj . The asymptotic stability of small amplitude solitary waves in the external potential 
has been studied in |Bou08[ IPSTO] . 

These models with self-interaction of local type have been receiving a lot of attention in the particle 
physics, as well as in the theory of Bose-Einstein condensates [HC091 |MJZ + 10 . The question of stability 



of solitary waves is of utmost importance: perturbations ensure that we only ever encounter stable configu- 
rations. Recent attempts at asymptotic stability |PS1Q[ IBC11] rely on the fundamental question of spectral 
stability: 

Consider the Ansatz ip(x,t) = (4>u>(%) + p( x > t))e~ lut , with 4> LLl (x)e~' lu ' t a solitary wave solution. 
Let dtp — A u p be the linearized equation on p. Does A u have eigenvalues in the right half-plane? 



In spite of a very clear picture of the spectral stability of nonlinear Schrodinger and Klein-Gordon 
equations [VK731 [Sha831 IWei85 , SS85, GSS87] and many attempts at the spectral stability in the context of 
the nonlinear Dirac (let us mention [AS83( IAS86| ISV86| IBou08| ICKMS10] ), in the latter case the question 
of spectral stability is still completely open. Our numerical results }BC09] show that in the ID Soler model 
(cubic nonlinearity) all solitary waves are spectrally stable. Let us mention the related results [CP06 , C hu07] . 

Our previous result [Comllj shows that the Vakhitov-Kolokolov criterion in the case of the nonlinear 
Dirac equation gives a less definite answer about the spectral stability than in the case of the nonlinear 
Schrodinger equation. All we know is that when d^Q(ui) = 0, with Q(u>) being the charge of the solitary 
wave (j) u e~ luJt , two eigenvalues collide at A = 0, but we do not know where these eigenvalues are located 
when d u Q(u>) 7^ 0. On the other hand, for the solutions to the nonlinear Schrodinger equation, the condition 
9ujQ{uj) > is enough to conclude that one positive and one negative eigenvalues move out of A = along 
the real axis; see |VK73l[GSS87l[CP03] . 

In the present paper, we show that if the Vakhitov-Kolokolov guarantees linear instability for the system 
obtained in the nonrelativistic limit, then the same result is true for the small amplitude solitary waves in the 
original relativistic system. We use our approach to show that the small amplitude solitary wave solutions 
to the nonlinear Dirac equation in ID with higher order nonlinearities are linearly unstable. 

According to our results, the spectrum of the linearization at small amplitude solitary waves in ID 
nonlinear Dirac equation with the cubic and quintic nonlinearities has no real eigenvalues; instead, one can 
prove the existence of two purely imaginary eigenvalues. To prove that these solitary waves are spectrally 
stable, one also needs to prove that there are no complex eigenvalues with Re A > 0; we already have partial 
results which we will publish elsewhere. 

Our approach is based on the idea that the family of real eigenvalues of the linearization of the nonlinear 
Dirac equation bifurcating from A = is a deformed family of eigenvalues of the linearization of the nonlinear 
Schrodinger equation. The model and the main results are described in Section [2] The necessary construc- 
tions in the context of the nonlinear Schrodinger equation are presented in Section [3l The asymptotics of 
solitary waves of the nonlinear Dirac equation and the linearization is covered in Section @] The statement 
of Theorem 12.11 for u> < m follows from Proposition 15.21 which we prove using the Rayleigh-Schrodinger 
perturbation theory. 

ACKNOWLEDGMENTS. The author is grateful to Gregory Berkolaiko, Nabile Boussaid, Maria Esteban, 
Dmitry Pelinovsky, Iosif Polterovich, Bjorn Sandstcdc, Walter Strauss, Boris Vainberg, and Michael Wein- 
stein for most helpful discussions. The author is grateful to Stephen Gustafson for the preprint |GG] with a 
similar result for the nonlinear Dirac equation in 3D. 

2 Main result 

We consider the nonlinear Dirac equation in one dimension, 

id t ip = D m ip - fWM)fal>, zeR, ^eC 2 , (2.1) 

where D m is the Dirac operator: 

D m = —iad x + rri/9, m > 0. 

Above, tp* being the Hermitian conjugate of ip. We assume that the nonlinearity /(s) is smooth, real-valued, 
and satisfies /(0) = 0. The Hermitian matrices a and j3 are chosen so that 

D 2 m = (- iad x + /3m) 2 = {-d 2 x + m 2 )I 2 , 

where I 2 is the 2x2 unit matrix. That is, a and /3 are to satisfy 

a 2 =I 2 , P 2 =h; a(3 + f3a = 0. (2.2) 

The generalized massive Gross-Neveu model, or, in the terminology of [CKMS10 , the scalar-scalar case with 
k G N, corresponds to the nonlinearity f(s) = s k . 



According to the Dirac-Pauli theorem (cf. |Dir28[ lvdW321 IPau36j and |Tha921 Lemma 2.25]), the par- 
ticular choice of a and /3 matrices does not matter. We choose 

For a large class of nonlinearities /(s), there are solitary wave solutions of the form 



ip(x,t) = 4>uj(x)e 



V(;U) 

u(-,oj) 



eH\R,C 2 ), \u\<m, (2.3) 



with v(x,lo) and u(x,uj) real- valued and with finite H 1 norms. For example, this is the case for the Soler 
model with the nonlinearity f(s) — s. For details, see Section @] 

Due to the U(l)-invariance, for solutions to (|2.ip the value of the charge functional 

Q(V)= [\i>(x,t)\ 2 dx 

Jr 

is formally conserved. For brevity, we also denote by Q(u>) the charge of the solitary wave 4> UJ (x)e^' lut : 

Q(u)= [ \M*)\ 2 dx. (2.4) 



We are interested in the spectral stability of linearization of (|2.1j) at a solitary wave solution (|2.3j) . 

Theorem 2.1. Let 

f{s)=as k +0{s k+1 ), a>0. 

If k > 3, then the solitary wave solutions <p u &~ liJt to (|2.1[) are linearly unstable for u € I, where I C Q 
is the largest interval with supp/ = m such that d u Q{uj) does not vanish on I. More precisely, let A^ be 
the linearization of the nonlinear Dirac equation at a solitary wave <j>u(x)e~ . Then there are eigenvalues 
±A W G o- p (A u ), with A w > for u> € I, with X u — 0(m — u>). 

See Figured! 

Remark 2.2. We only need to prove the linear instability for lu < m. Then, by |Comllj . the positive and 
negative eigenvalues remain trapped on the real axis, not being able to collide at A = for u) E I = (u>o, rn) 
as long as d u Q{uj) does not vanish on / (this nonvanishing is a sufficient condition for the dimension of the 
generalized null space to remain equal to two). These eigenvalues can not leave into the complex plane, 
either, since they are simple, while the spectrum of the operator is symmetric with respect to the real and 
imaginary axes. 

Remark 2.3. We do not know what happens at u> = inf /: as lu drops below inf /, it could be that either the 
pair of real eigenvalues, having collided at A = 0, turn into a pair of purely imaginary eigenvalues (linear 
instability disappears), or instead two purely imaginary eigenvalues, having met at A = 0, turn into the 
second pair of real eigenvalues (linear instability persists). 

Remark 2.4. Theorcm l2.1l is in the formal agreement with the Vakhitov-Kolokolov stability criterion |VK73j . 
since for ui < m one has Q'(lo) < for k = 1 and Q'(w) > for k > 3. Let us mention that the sign of the 
instability criterion, Q'(ui) > differs from [VK73] because of their writing the solitary waves in the form 

tp{x)e +lut . 

Remark 2.5. We expect that in the ID case with k = 1 the small amplitude solitary waves are spectrally 
stable; we will prove this elsewhere. 

We also expect that in the ID case with k = 2 ("quintic nonlinearity", / (ip* f3i/j) f3ip = 0(\ip\ 5 )) the small 
solitary wave solutions of the nonlinear Dirac equation in ID are spectrally stable. For the corresponding 
nonlinear Schrodinger equation (quintic nonlinearity in ID), the charge is constant, thus the zero eigenvalue 
of a linearized operator is always of higher algebraic multiplicity. For the Dirac equation, this degeneracy 
is "resolved": the charge is now a decaying function for uj < m (with nonzero limit as u> — > m), suggesting 
that there are two purely imaginary eigenvalues ±A W in the spectrum of A u , with A w = o(m — to), but no 
eigenvalues with nonzero real part. 



i(m + lo) 
i(m — lo) 



K 



Figure 1: Main result: The point spectrum of the linearization of the nonlinear Dirac equation with the nonlinearity 
f(s) — s k + 0(k k+1 ), k > 3, at a solitary wave with lo < m contains two nonzero real eigenvalues, ±A^, with 
A„ = 0(ra — lo). See Theorem 12.11 Also plotted on this picture is the essential spectrum, with the edges at 
A = ±i(m — lo) and with the embedded threshold points (branch points of the dispersive relation) at A = ±i(m + ui) . 

Remark 2.6. The same approach can be used to show instability of small amplitude solitary wave solutions 
to the nonlinear Dirac equation in n > 3, with f(s) — as k + 0(s k+1 ), a > 0, with any k £ N. These results 
in 3D have been independently obtained in |GGj . Let us notice that in the 3D case for the cubic nonlinearity 
f(s) = s (this is the original Soler model from |Sol70j ). based on the numerical evidence from |Sol70, AS83 , 
one expects that the charge Q{ui) has a local minimum at u> = 0.936m, suggesting that the solitary waves 
with 0.936m < u> < m are linearly unstable, but then at u> = 0.936m the real eigenvalues collide at A = 0, 
and there are no nonzero real eigenvalues in the spectrum for ui < 0.936m. 

Remark 2.7. We can not rule out the possibility that the eigenvalues with nonzero real part could bifurcate 
directly from the imaginary axis into the complex plane. Such a mechanism is absent for the nonlinear 
Schrodinger equation linearized at a solitary wave, for which the point eigenvalues always remain on the 
real or imaginary axes. At present, though, we do not have examples of such bifurcations in the context of 
nonlinear Dirac equation linearized at a solitary wave. 

Remark 2.8. The existence of solitary waves stated in the theorem follows from |CV86] . We will reproduce 
their argument in order to have the asymptotics of profiles of solitary waves for lo near m. 

3 Nonlinear Schrodinger equation 

We are going to use the fact that the nonlinear Dirac equation in the nonrelativistic limit coincides with the 
nonlinear Schrodinger equation, 



id t il> = --dli> + rmjj- J{\iP\ 2 )iii ijj(x,t)£C, x£ 



(3-1) 



We will assume that 



f(s) = s k , k£ N. 



Solitary waves 

The solitary wave solutions 



ip{x,t) = 4>u,(x)e~ 



^ G H\R), 



exist for ui G (—00, to). These solitary waves and their asymptotics are well-known, to the extent that we do 
not know which of the multiple references would be most appropriate. Often we will not indicate explicitly 
the dependence of the amplitude of the solitary wave on ui, writing (j> instead of <p u . This amplitude satisfies 
the equation 

d 2 x 4>{x) = 2(m-u)(f)-2(f) 2k+1 , lei, 

which could be integrated to the relation 



d x (t> = -<f>\ 2 (to - u) 



20 



2 k 



jfc+1 



Introducing 3£(x, ui) = <^,(x), we get 



d x % = -23£\2{m-u)) - 



25£ k 
fc + 1' 



We will perform the scaling in terms of 
Then, as can be seen from (|3.2|) . 



y 2(m — ui). 



&(x,uj) = E 2/k U(ex), 
with U(y) positive spherically symmetric (even) solution to the equation 



(3.2) 

(3.3) 

(3.4) 



d y U(y) = -2EM1 - ^j^M, lirn U(y) = 0. 

Such a solution exists and is unique; it is explicitly given by 

A- + 1 \ i/fc 



U(y) 



2 cosh ky 



(3.5) 



Linearization at a solitary wave 

To derive the linearization of the nonlinear Schrodinger equation (|3.1[) at a solitary wave ij)(x, t) — W (x)e~ 
we use the Ansatz 

V>0M) = {Mx)+P(x,t))e~ iut , p{x,t) G C, 

and arrive at the linearized equation 



x G 



d t p = )L(uj)p, p(x,t) = 



Rep(x,t) 
Jmp(x,t) 



with 



1 
-1 



L_ 

L+ 



L_(w) = --81 +m-u- M), L + (uj) = L_- 2f'{tf,)<&. 



(3.6) 



(3.7) 



(3.8) 



In the case of the nonlinearity f(s) — s , taking into account the explicit form of 3£{x) given by p.4[) , one 
obtains 



L-{w) = --B% + m-u-?' 



fc + 1 



1„, ,(2fe + l)(fc+l) . . 

i + H = --^ + TO- w -e 2 ^- 4 , -i . (3.9) 



2 cosh (efcx) 2 2 cosh (sfcx) 

Lemma 3.1 (Vakhitov-Kolokolov stability criterion). For the linearization (J3.6I) a< a particular solitary 



wave <f) ul (x)e lujt , there are real nonzero eigenvalues ±A G o~d(jL), A > 0, if and only if ■£j\\<t>w\\ 2 ^2 > at 



this value ofui. 



Let us consider the spectrum of the linearized equation (13.61) for the nonlinearity f(s) = s k . Using 
one derives for the corresponding charge: 



Q{u) 



(x)dx= / 3E{x,w)dx = e* l C = (2(m- w))* =C, w < m, 



(3.10) 



where C = J* R [/(y) dy > 0. We see from (J3.10I) that one has Q'{u) < for fc = 1, Q'{ui) = for k = 2, and 
Q'(u>) > for fc > 3. 

Lemma 3.2. Let /(s) = s fc ; fc € N. If k = 1 or fc = 2, i/jen a p (jL) C iR. //fc > 3, i/ien c7 p (jL) 9 {±e 2 A}, 
/or some A > and e = y/2(m — u). 



Proof. In the case fc = 1, since Q'(uj) < by p.lOj) . Lemma T3.ll guarantees that there are no nonzero real 
eigenvalues: <j p {)L{uj)) n (K\0) = 0. In the case k = 2, Q'{ui) = 0, and the eigenvalue A = is always of 
increased algebraic multiplicity 6 (generically, the algebraic multiplicity of A = is 4, jumping to 6 when 
Q'(ui) = 0). In the case fc > 3, since Q'(w) > by p.lOj) . Lemma [3TT1 states that there are two real eigenvalues 
±A U € (T p (jL(u)), with A w > 0. The rescaling y — ex shows that A w = e 2 A, with some A > 0. □ 

For fc € N, consider the Schrodinger operators 



L_ 

L+ 



L_ 



1 l__fc + l_ 

2^ + 2 2cosh 2 fcy' 



L + = -\ d l 



1 (2fc + l)(fc+l) 

2 2 cosh 2 fcy 



(3.11) 



Let 4>(y) = 



cosh 1 / fc ky ' 



then 



L-cf) = 0, L+d x <$> = 0. 



Proposition 3.3. If k > 3, £/ie discrete spectrum of the operator )L 

nonzero eigenvalues ±A, A > 0. If k <2, o~d(jL) C zK. 

Proof. We start with constructing several important relations. By (|3.12l) . 

fc(fc + 1 



L_ 

-L+ o 



(3.12) 
contains two real 



L+4) = L_4) + (L + -L_)4) 



cosh ky 



(3.13) 



For fj, > 0, let cb M (y) = cb(^/) and L M = -±<9 2 + ^i 2 ( \ - j^pr^J ' then ^n^n = °- Taking the derivative 
of this relation with respect to fi and evaluating it at \x = 1, we get: 



L,0 



fc + 1 
cosh ky 



*, 



where 9 



9/x 



1 y sinh fcy 

— iiji 

i= iCOsh fc /ifcy cosh k ky 



By (l3~T3l) and (|3Tl]l . 



L+(-9 - -4)) = 4). 



(3.14) 



(3.15) 



Now we follow the argument from [VK73] . 



Lemma 3.4. Let fc > 3. T/ie minimum of (r, L + r) under the constraints (r,r) = 1 and (4>,r) — is 
negative. 

Proof. The vector r corresponding to the minimum of (r, L+r) under the constraints (r, r) = 1 and (<J), r) = 
satisfies the equation L + r = ar + /3cb, where a, j3 G K are Lagrange multipliers. Since (r, L+r) = (r, or + 
/3cf>) = a, we need to know the sign of a. Denote 



f(z) = (4>,(L+-z)- 1 4>), ze P {L + ). 



We note that f(z) has a removable singularity at z — since 4> £ L 2 (R), while the restruction of L + onto 
the space of spherically symmetric (that is, even) functions, L + : L 2 (R) — > L 2 (R), has a bounded inverse, 
(being even) Using (|3.15|l . we have: 



/(0) = (4>,i; 1 4>> = (4),(- 

/ z sinh z dz \ dz 

v cosh fc z cosh fe z 



:*)) 



-zd 



1 



?/ sinh fcy 



1 



dy 



cosh fc ky V cosh k ky /ccosh fc fcy 

1 \ dz \ , 1 1 \ /" <^ z 

2 J I J ~ U — TJ I I 

2 cosh fc z y A: cosh k z ' L K Jr cosh fc z 



>0, 



where we took into account that k > 3. Since /(0) > and lim /(z) = — oo (it is — oo since f'(z) > 0), 

z->A 4- 



where Aq < is the smallest eigenvalue of L + , there is /x £ (Aq, 0) such that /(/x) = 0. 



D 



Since L_ is positive-definite and r in Lemma 13.41 is orthogonal to kerL_ (which is spanned by it)), we 



r — 1/2 

may define it = L_ r; then 



{R,L l J. 2 L + l]l 2 R) = (r,L+r) < 0. 

Therefore, ad(L-L+) — <Jd{L_ L + LJ ) contains a negative eigenvalue —A 2 , where A > 0. Let £ be a 
corresponding eigenvector, so that L-L + £ = — A 2 £. Then 



£_ 

-L+ 



' AC ' 



±A 



AC ' 

T^4-£ 



hence ± A £ aa(jL). 



D 



Remark 3.5. Since the spectrum of the linearization at zero solitary wave is purely imaginary, one has 
lim X u — 0. 



Remark 3.6. Comparing L (see (|3.7[) . (|3.9[1 ) to the operator L introduced in p. tip . 

fc + 1 



L_ 

L+ 

one concludes that 



2 y 2 2cosh 2 fc2/ 



^4 



1 1 (2 fe + l)( fe + l) 

2 y 2 2cosh 2 fczj V ' 



a(L-)=e 2 a(L^), a(L+) = e 2 a(L+), a(jL) = e 2 a()L). 



Therefore, A w = e 2 A, where A is the positive eigenvalue of )L. 



4 Nonlinear Dirac equation 

In this section, we will use the notation 



In terms of the components, ip 



^2 



g(s) = m$-f(s). 
■01, tp2 S C, we write the nonlinear Dirac equation (12.11) as a system 



zcWi = a E </>2 + s(IV'il 2 -IV'2| 2 >/'i 



id t i>2 



-SxVl" 5(1^1 -1^2 1 2 )^2 



(4.1) 



Solitary waves 

Definition 4.1. The solitary waves are solutions to (|2.1j) of the form 

S = {ip(x,t) = (f> LJ (x)e~ luJt ; <$>„ G H\R,C 2 ), u G R}. 

We start by demonstrating the existence of solitary wave solutions and exploring their properties. The 
following result follows from |CV86j : we follow our article [BC09 J. 

Lemma 4.2. Let G be the antiderivative of g such that G(0) = 0. Assume that there is loq < m such that 
for given to G (ojq, to) there exists r u > such that 

ur u = G{r u ), uj ^ g(r u ), and us < G{s) for s G (0, r u ). (4.2) 

Then there is a solitary wave solution ip(x,t) = <f> u (x)e~ to (|2.1|) . where 



<t>u{x) 



v(x,Uj) 

u(x,u) 



(4.3) 



(4.6) 



with both v and u real-valued, belonging to H (R) as functions of x, v being even and u odd. 
More precisely, for x G R anrf w G (o; ,to), Zei us define 3£{x,w) and W(x,w) by 

% = v 2 - u 2 , <¥ = to. (4.4) 

T/ien ^(xjw) is a unique positive symmetric solution to 

d 2 x X = -d x (-2G(3£) 2 + 2u 2 S£ 2 ), lim ar(ar,w)=0, (4.5) 

x— >±oc 

and ^(x,w) = — J^c? xt $K"(x, w). TTiis solution satisfies &(0,w) = r u . 

Proof. From (|4.1[) . we obtain: 

UW = c^u + .<?(H 2 — M 2 )f, 
uju = —d x v — .g(|u| 2 — |u| 2 )u. 

Assuming that both v and u are real- valued (this will be justified once we found real- valued v and u), we 
can rewrite (|4.6p as the following Hamiltonian system: 

d x u = u>v- g{v 2 ~u 2 )v = d v h(v,u), , , 

—d x v = uju + g(v 2 — u 2 )u~duh(v,u), 

where the Hamiltonian h(v,u) is given by 

h(v,u) = ~(v 2 + u 2 )-±G(v 2 -u 2 ). (4.8) 

The solitary wave with a particular w G (wo, m) corresponds to a trajectory of this Hamiltonian system such 
that 

lim v(x,lo) = lim u(x,uj) = 0, 

x— »±oo x— v±oo 

hence lim a ._>.-|- 0o i£T = 0. Since G(s) satisfies G(0) = 0, we conclude that h{v(x),u(x)) = 0, which leads to 

uj{v 2 +u 2 ) = G{v 2 -u 2 ). (4.9) 

We conclude from (J4.9I) that solitary waves may only correspond to |w| < m, w ^ 0. 
The functions 3£{x,u)) and £^(x,w) introduced in (|4.4|) are to solve 

d x W = -{v 2 + u 2 )g{&) + u % = -±G(3£)g(3£) + wJT, [ UJ 

and to have the asymptotic behavior limui^oo JT(x) — 0, limui^oo W{x) — 0. In the second equation in 
(|4.10[) . we used the relation (14.91) . The system (|4.10p can be written as the following equation on 2£\ 

d 2 x X = -%-(-2G( St) 2 + 2u?3t 2 ) = 4(G(&)g(&) - w 2 JT). (4.11) 

This equation describes a particle in the potential — 2G(s) 2 + 2w 2 s 2 . The condition (J4.2I) is needed for the 
existence of the turning point of the zero energy trajectory in this potential, at $ = r u . □ 



Solitary waves in the nonrelativistic limit 

From now on, without loss of generality, we will assume that the nonlinearity / £ C°° (R) is such that 

f(s) = s k + o( S k+l ), k£N, 

and consider k fixed throughout the rest of the article (this is to avoid writing subscripts). By Lemma l4.2| 
there is ujq < m such that there are solitary wave solutions for u) £ (u>o,m). 

We are going to determine the asymptotics of solitary waves in the nonrelativistic limit u — > in. 

We consider the function g{s) := m — f(s) and its antiderivative, 

g(a) = to- s k + 0{s k+1 ), G(s) =ms- + 0(s k+2 ). 

k + 1 

Then the expression in right-hand side of (J4.11I) takes the following form: 

G(s)g(s) - lo 2 s = (ms - ^- + 0( S k+2 )) (to - s k + 0(s k+2 )) - u 2 s = (m 2 - u 2 )s - 1 l±l s k +^ + 0(s k+2 ). 

Let u> £ (loq, m), so that there is a solitary wave 4>cj(x)e~ lut with this particular value of uj. Denote 



e = \/m 2 ~uj 2 . (4.12) 

Let 3^(x,(jj) be a positive symmetric solution to (|4.11|) from Lemma T4.2( we rewrite (|4.11j) as 

-i^ t r-^^s" fc+1 + e 2 jr = o( I r fc+2 ). (4.13) 

4 x k + 1 y ' y ' 

Let U(y) be a positive symmetric solution to 

- - A dp - |^V +1 + U = 0, km U{y) - 0. (4.14) 

4 J k + 1 y^oo 

Such a solution U(y) exists and is unique, and is explicitly given by (|3.5p . 

Lemma 4.3. For ui £ (uj ,m), there are the relations 

3T(x,uj) = eiu{ex)+0(eT;), (4.15) 

v(x,u) = ei(U(ex))7 +0(e%), u{x, uj) = 0(e 1+ i). (4.16) 

Remark 4.4. Similar asymptotics are establisked in .Gua08] for tke nonlinear Dirac equation in 3D. Tkese 
asymptotics are used in }GG| for tke proof of linear instability of solitary waves witk u> < m to tke nonlinear 
Dirac equation witk cubic nonlinearity. 

Proof. Since U{y) > for y £ K, one concludes from (|4.14p tkat A — is tke lowest eigenvalue of tke 
operator 

H - = --A-r ± ^ uk + l - 

4 v k + 1 
Taking tke derivative of (|4.14l) . we find: 

H + d y U := -jd%U - (k + 2)U k d y U + d y U = 0. (4.17) 

Tke operator H + in tkis relation is tke same as in (I3.11[) . Witk d y U kaving one node, we conclude tkat 
A = is tke second lowest eigenvalue of H + . 



We define the function X(y, e) by the relation 2£{x, w) = e 2 / k X{ex, e), where e and lj are related by 
flUTg) . Then (J4TT3]) takes the form 

- i^ 2 X + X - ^±^ fc+1 = e- 2 -*0( JT fc + 2 ) = £ 2 / fc O(X fc + 2 ). (4.18) 

4 y fc + 1 



Subtracting (J4.14I) from (J4.18I) . we find: 



_I 3 2 (X - 17) + (X- 17) - Y^(X k+1 - U k+1 ) = e 2 ^ k O(X k + 2 ), 

which we rewrite as 

-\d 2 {X -U) + (X-U)-(k + 2)U k (X -U) = e 2/k O(X k+2 ) + 0((X - U) 2 ). 
4 y 

Thus, denoting Z(y,e) = X(y,e) — U(y), one has 

Z= (H + \ L2 y 1 (e 2 ' k O{U k+2 )+€ 2 ' k O{Z) + 0{Z 2 ) S j. (4.19) 

Above, H+ 1 2 : L 2 (M.) — > L 2 (R) is the restriction of H+ to the space of spherically symmetric (even) functions; 

its inverse is bounded (from L 2 to H}) since by (|4.17[) we know that kerH + is spanned by d y U , which is 
skew-symmetric (let us mention that A = is a simple eigenvalue of H + , which is straightforward in one 
dimensional case). The right-hand side of (J4.19I) is well-defined since both U and Z = X — U arc spherically 
symmetric. 

For given e > 0, let B e i/u = {Z £ # X (R); Z(x) = Z(-x), \\Z\\ H i < e 1 ^} C ff*(R) be the space of 
symmetric functions with ||Z||#-i < e l ' k . Consider the map 

S e i A -> H^QBL), Z4 (H+l^y 1 (e 2/k O(U k+2 ) + e 2/k O(Z) + 0(Z 2 )Y (4.20) 



We can choose uii £ (uiq, to) so close to m, that for ui £ (wi, to) the value of e = y/m 2 — to 2 is small enough 
for the map (J4.20I) to be an endomorphism and a contraction on B e i/h] by the contraction mapping theorem, 
there is a unique stationary point Z £ B e x/k, and, by f)4.19[) . Z satisfies ||Z||#-i = 0(e 2 ^ k ). It follows that 

\\3r(x,u) - e 2 ' k U{tx)\\ L oo = e 2 l k \\Z\\ L ~ m = 0(e 4 / fc ). 

Now the asymptotics of ^(a;,w) could be determined from (J4.10I) and then v{x,uj), u(x,ui) are determined 
by (gU). ' D 

5 Linear instability of small solitary waves 
Linearization at a solitary wave 



To derive the linearization of equation (|2.ip at a solitary wave (|2.3[) . we consider the solution in the form of 
the Ansatz 



ij{x,t) = ((f> u (x)+p(x,t))e- iut , u (a;) 



v{x,uj) 
u(x,uj) 



Gff 1 (R,C 2 ), (5.1) 



where p(x,t) £ C 2 . Note that v(x,u>) and u(x,w) are real-valued by Lemma 14.21 Then, by (|2.1[) . the 
linearized equation on p takes the following form: 

ip = D m p -up- f(<f>Z0<l> u )Pp - (fcPp + P*P<t>u)f'((f>ZP(f>u,)l3<t> u , D m = -iad x + to/3. (5.2) 

We note that the above equation is R-linear but not C-linear, due to the presence of the p* term. We denote 
, with Rj, Sj real- valued. Since v and u are real- valued, 



P 



R 2 + iS 2 



4>lPp + p*H^ = 2 Re(Op) = 2(wfli - uR 2 ), 
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we rewrite equation (J5.2I) in the following form in terms of R 



Ri 
R 2 



di 



where 

and J corresponds to 1/i: 



D m -w-fP \R 

D m -u-fp S 



f = f(cf>ZP4> u ), f = f 



h 
-h o 



c _ S\ 
. 2 . 

2Re(C/W 



(5.3) 



where /2 is the 2x2 unit matrix. We define L±(w) and L(w) by the following: 



L_(w) = 



m — u — f d x 

-d x -m — u + f 



L+M 



',.2 



to — u; — / — 2/ t; 

-9a: + 2/'vU 



& + 2/'im 



-m — u> + f — 2/ m 



',,2 



(5.4) 



L(w) 



L+(w) 
L_(w) 



Let us remind the reader that both u and it in (15.41) depend on u>. Then equation (15.3[) which describes the 
linearization at the solitary wave 4> u e~ lu)t takes the form 



9 



JL(«) 



I? 







L_M 




R 


5 




-L+(«) 







S 



wave u e 



Lemma 5.1. For any nonlinearity /(s) in (|2.ip . t/ie spectrum of the linearization at a solitary 
has the following properties: 

1. a ess (JL(w)) = iR\(-i(m - u),i(m - w)); 

& o-«j(JL(w)) D {±2w«; 0}. 

See [Comllj . 

Unstable eigenvalue of JL for u < m 

Proposition 5.2. Let k > 3. There is ujq < m such that for u G (wo, to.) i/iere are iwo families of eigenvalues 

±A W G o- p (JL(w)), A w > 0, A w = 0(e 2 ). 



Proof. In the explicit form, the relation 

-A 



L_ 

-L+ 



can be written as follows: 



g — ijj d x 

—A —d x —g — u 

m + uj + f + 2/V -d x - 2f'vu -A 

d x -2f'vu m + uj- f + 2f'u 2 -A 



'Ri 

Si 
S 2 



= 0. 



(5.5) 



Dividing the first and the third rows by e 2 , the second and the fourth rows by e, and substituting y = ex, 
i?2 = eR.2, S2 = eS2, we get 



A 

~77 





-m+uj+f+2f'v 2 



d y -\2f'vu ]{m + co- f + 2f'u 2 ) 
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lg _ 2£i 

- ^7/ -2 



e 2 


^v 




«1 


-9„ 


i(_m-« + /) 




e-R 2 


A 
~1? 







5i 





_ A 




e<S Y 2 



(5.6) 
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which simplifies to 



A 

~7? 





-A 

-d -2£ 



m—ui—f 

3 



',,2 





-m+uj+f+2f'v 2 



d y - ±2fvu m + uj- f + 2f'u 



Let A = lim -4. We introduce the matrices 



k-A = 



By 

—m — u + f 




~R{ 
4 







Si 


-A 




IS2] 



(5.7) 



-A 


\ - u k 


0„ 





-dy 


-2 


-\ + {2k + l)U k 


~dy -A 





By 


2 






Ki =diag[l, 0,1,0], K 2 = diag[0, 1,0, 1]. (5.. 



Above, 



cf. 



U(y) = 



A: + 1 \ Vfc 



• 2 cosh ky 

Noting that one has 3£(x) = e*U(ex) + O(e^) by Lemma T4.3[ we can write (|5.7p in the form 

A 



(5.9) 



Aat? = (-j - a)K!T? + AK 2 ?7 + WV7, 



(5.10) 



where r\ = 



Ri 

Rz 

Si 

s 2 



G C 4 and 



W(y,e)=Ai. 



A 

"77 



m—OJ — f 
71 







-m+uj+f+2f'v 2 


dy 


2f'vu 

6 



-X 

_ a _ 2f'vu 

U V £ 



-d 



dy 

-m — uj + / 



-A 



-(A_ A )Ki-AK 2 (5.11) 



m + u - / + 2/V 

is the zero order differential operator with L°° coefficients. 
Lemma 5.3. ||W(-,e)|U» (Ri c*_ > c*) < 0{e 2 ' k ). 
Proof. By Lemma |4.3[ one has 

v(x,u) = e*(J7(er))* + 0(e*), IK,*- = 0(e 1+ i). 

Then 

/( V 2 - u 2 ) = v 2k + 0{e 2+ i) = e 2 U{ex) k + 0(e 2+ *), 

f(v 2 - u 2 ) = kv 2k - 2 + 0{e 2 ) = (9(e 2 -*), 
f'{v 2 - u 2 )v 2 = e 2 kU{ex) k + 0(e 2+ %), 
f'(v 2 -u 2 )vu = 0{e 3 ). 
Now the proof follows from the definition of Aa and Ki , K 2 in (|5.8[) . 
Lemma 5.4. dimker Aa > if and only if A is an eigenvalue of the operator 



□ 



JH 



\dl - \ + {2k + l)U k 



-Id 2 + i - t/ A 

2 c/ y -r 2 o 





When k > 3, t/iere is A > smc/i iftai ±A e er rf (jH). 
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Proof. The first statement of the Lemma follows from the structure of the second and the fourth rows of 
Aa- The second statement follows from Proposition 13.31 □ 

By Lemma [5.41 if k > 3, then there is A > such that ±A G Od QH). Now we will use the Rayleigh- 
Schrodinger perturbation theory to show that there are A G <7</,(JL) with A ra ±e 2 A. Let $a G kcrAA 
with ||$a||l 2 = lj an d let Pa : L 2 (R, C 4 ) — > i7°°(R, C 4 ) be the spectral projector onto the corresponding 
eigenspace. Coupling (|5.10|) with J$a and projecting (|5.10|) onto Range (1 — Pa), one has: 



= < J$ A , A AV ) = ( A _ A ) ( J<T> \ . K i. // ) - , A ( J<T> A •• K 2 ,/} + .; J<T> v . M - //} . 



(5.12) 



(1 - P A )r? = A^ x (l - Pa) ((^ - a)Ki + AK 2 



W h? 



We will find 77 in the form 77 = $a + C with £ G Range (1 — Pa). Denote // 



(5.13) 
A. Then the relations 



(|5.12p and (J5.13I) can be written as 

= /i(J$ A , Ki* A ) + ^(J$a, KxC) + e 2 (A + m)<J*a, K 2 ($a + C)> + (J*A, VF($ A + C)>, (5-14) 

C = AX : (1 - Pa)(/*Ki + e 2 (A + /i)K 2 + w) ($ A + 0- (5-15) 



Let us note that (J$a, Ki$a) 7^ 0. The equations (J5.14I) . (I5.15J) could be written as /i = M{n, (), ( — Z(/i, (), 
with functions M:lx L 2 (R, C 4 )^M,2:Mx L 2 (R, C 4 ) -> L 2 (R, C 4 ) given by 



M(/*,0 = - 



1 



(J$ A ,Ki$/ 



/i(J$ A , KiC) + e 2 (A + /i)(J$ A , K 2 ($ A + 0) + (J$A, ^($ A + 0) 



Z{ji, C) = A^l - Pa) UKi + e 2 (A + M )K 2 + W) (4 A + C) 



Pick T > 1 such that 



T> 2||AX 1 (l -Pa)Ki* a |U». 



(5.16) 

(5.17) 
(5.18) 



(5.19) 



Lemma 5.5. Consider R x L 2 (M,C 4 ) endowed with the metric 

IKM.OIIr =r|/<| + ||CIIl2(e,c*)- 

There is Wi £ (woj w) suc/i t/iai /or w € (wi, m) t/te map 

MxZ: Ix J L 2 (R,C 4 )^Rxi 2 (M,C 4 ), (fi,Q >-> (M(^,C), £(m,C))> 

restricted onto the set 

B el/k ={(/i,C) GKxL 2 (R,C 4 ); ||(/J,C)l|r < e 1/fc } C R x L 2 (K,C 4 ) 

is an endomorphism and a contraction with respect to \\ ■ ||p. 

Proof. The proof follows from the definitions (|5.16j) . (|5.17j) . We needed to introduce the factor T into the 
definition of the metric to make sure that the contribution of the term A^ 1 (l — Pa^Ki^a from (|5.17|) 
into ||(Af(/j, Q,Z(n, C))||r is bounded by ^|/i| < |||(/J, C)l|r- The contribution of all other terms is bounded 
by j||(M>C)l|r + 0{e 2 / k ), placing {M{ix, Q,Z(/j,,Q) inside B e i/k whenever (//, Q G B e i/k, as long as wi is 
sufficiently close to m, so that e = \fmF~~u? is sufficiently small. The contribution 0(e 2 / fc ) comes from the 
terms with VT$a, due to the bound \\W(-, £)||l°°(r,c 4 ^c 4 ) = 0(e 2 / fe ) established in Lemma [5731 

The contribution of the term A^ : (l — Pa)m^i < I ) a into the norm \\Z(fi, £) — Z(/j>' ,(')\\ L 2 is bounded by 
±r|/i-//| < i||(/j-/i',C-C')llr- The contribution of all other terms from ([5TTB| . ([53?) into ||(M(/x,C)- 
M(//,£'), Z(n, () — Z (fi' , ( r j)\\r could be made smaller than |||(/U — /J,', C — Ollr by choosing wi sufficiently 
close to m, so that e = Vw 2 — w 2 is sufficiently small. It follows that (M x Z)\ g is a contraction in the 
metric II ■ lip. □ 
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According to Lemma 15.51 by the contraction mapping theorem, for u> £ (wi,m), the map (J5.19I) has a 
fixed point 

(fx (oj), Co(«)) G B 6 v* C R x L 2 (M, C 4 ). 

Thus, we have 

±e 2 (A + /^(w)) e a p (JL(w)), w G (wi, m), 

with r|/xo(w)| < e 1 /' , finishing the proof of the proposition. □ 

By Remark 12.21 Proposition 15.21 finishes the proof of Theorem 12. II 
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